Possible magnetic-field-induced voltage and thermopower in diluted magnetic semiconductors 
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In diluted magnetic semiconductors, the carrier concentration and the magnetization of local moments are 
strongly coupled, since the magnetic interaction is mediated by the carriers. It is predicted that this coupling 
leads to an electric polarization due to an applied magnetic-field gradient and to the appearance of a magnetic- 
field-dependent voltage. An expression for this voltage is derived within Landau theory and its magnitude is 
estimated for (Ga,Mn)As. Furthermore, a large contribution to the thermopower based on the same mechanism 
is predicted. The role of fluctuations is also discussed. These predictions hold both if the magnetization is 
uniform and if it shows stripe-like modulations, which are possible at lower temperatures. 
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I. INTRODUCTION 

Among magnetically ordered materials, diluted magnetic 
semiconductors 1 ' 2 ' 3 (DMS) are unique in the possibility of 
tuning the interaction between local moments in situ. Sev- 
eral groups 4,5 - 6 have demonstrated that the magnetization and 
Curie temperature T c in DMS can be changed by applying a 
gate voltage in a field-effect-transistor geometry. The gate po- 
tential affects the concentration of carriers, usually holes, in 
the DMS sample, which leads to a change of the magnetic 
interactions, since these are mainly mediated by the carriers. 
A strong dependence of the magnetization on carrier concen- 
tration is also seen in experiments on series of samples with 
different dopand concentrationsii^ 

On a microscopic level, each local moment magnetically 
polarizes the carriers in its vicinity due to the exchange inter- 
action between carriers and localized (usually d-shell) elec- 
trons of impurity ions. The other local moments see these 
carrier spin polarizations, leading to an effective magnetic in- 
teraction between local moments. This is the essence of the 
Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction^ It is 
here complicated by band-structure effects and spin-orbit cou- 
pling in p-type DMSpiii 2 *^ 4 ^ screening due to disorder, 15 
the small Fermi energy compared to the effective Zeeman 
splitting, 16 and competing short-range interactions. The 
RKKY-type interaction depends on carrier concentration es- 
sentially through the density of states at the Fermi energy. 
This dependence can be expanded to linear order in the con- 
centration, in agreement with experimental observations. 

The present author has recently suggested that this coupling 
may stabilize an equilibrium state with periodic, but strongly 
anharmonic modulations of carrier concentration and magne- 
tization below a certain temperature T* < T c pii The tem- 
perature T* strongly depends on rj = dT c /dn, the change 
of T c with carrier concentration. In Landau theory we find 
T c — T* oc 1/rj 2 . While some experimental results 8 are con- 
sistent with a small value of T c — T* ~ 10 K, which would 
make the effect observable, it is not clear whether this value 
is typical for DMS. In the present paper we put the main em- 
phasis on the uniform phase. 

In particular, we explore the reverse effect, compared to the 
gate control of magnetism: Can a change in magnetization 
lead to a change in carrier concentration? Or, if the sample is 



electrically isolated, can it lead to the appearance of a voltage? 
Microscopically the dependence of the effective spin interac- 
tion on the carrier concentration means that there is an ad- 
ditional contribution to the carrier energy depending on spin 
orientations. In particular, if the impurity spins are partially 
aligned, there is an attractive potential for carriers. This ef- 
fect is discussed and its magnitude is estimated in Sec. UJ 
In addition, the magnetization of course also changes with 
temperature, which should generate a contribution to the ther- 
mopower. This possibility is discussed in Sec.|n] The paper 
is summarized in Sec. lIVI 



II. MAGNETIC-FIELD-INDUCED VOLTAGE 

Let us consider the setup sketched in the inset in Fig.^ A 
metallic DMS sample is placed in an inhomogeneous mag- 
netic field and the voltage between opposite surfaces is mea- 
sured. We assume the temperature to be close to T c so that 
the magnetization is far from saturation. The magnetization 
is larger in the region with strong magnetic field. This region 
thus attracts carriers due to the mechanism discussed above. 
An electric current will only flow for a short time until the 
accumulated charge produces an electric field that prevents 
further charge accumulation. For positive (negative) carrier 
charge the electric potential in the strong-field region will be 
higher (lower) than in the weak-field region. The voltage mea- 
sured between the two regions divided by the magnetic field 
difference will thus have the same sign as the carrier charge. 

Following Ref. 17 we start from a Landau theory for the 
coupled magnetization and carrier density, defined by the 
Hamiltonian H = H m + Hs n with 
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The first part is the usual Landau functional for a Heisenberg 
ferromagnet in a magnetic field B directed along the z axis. 
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Summation over i = x, y, z is implied. The second part con- 
tains the screened Coulomb interaction and an applied electric 
potential. The carriers are assumed to have charge q = ±e. 
We impose the constraint J d 3 r 5n(r) — 0. The two parts are 
coupled through the coefficient a in H m , which depends on 
temperature and vanishes at the Curie temperature. We write, 
to leading order, a — a' (T — T c — r/ Sn), where T c is now the 
Curie temperature for Sn = 0. 

The equilibrium states can be found from the Euler equa- 
tions SH/5m = and 5H/S(Sn) = 0. 17 We find a param- 
agnetic phase for T > T c , a uniform ferromagnetic phase for 
T* < T < T c , and a phase with periodic modulation of m 
and Sn for T < T*, where T* = T c - e 2 f3j/e Q ea' 3 n 2 . By 
rescaling of length, magnetization, and energy it is possible to 
reduce the number of parameters. For example, the magneti- 
zation can be written in the scaling form 
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FIG. 1: Linear response of the voltage for a weakly nonuniform 
magnetic field, dV/dB oc m z dm z /dB, as a function of reduced 
temperature for various magnetic fields, assuming positive carrier 
charge. Scaling factors and reduced temperature are defined in 
Eqs. 0-^6). Inset: Sketch of DMS sample in a nonuniform field. 



where 




(4) 



are the natural units of magnetization and magnetic field, re- 
spectively, 



t = 



_ a'(T-T c ): 
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is a dimensionless reduced temperature, and 



fJ 



VJ3r t 
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(6) 



is a dimensionless measure of the coupling r/ between carrier 
concentration and magnetization. M is a dimensionless vector 
scaling function. 



A. Uniform phase 

In the uniform ferromagnetic and paramagnetic phases, we 
can drop the gradient term in H m . We obtain the energy den- 
sity h = h m + hsn with 



hs-n = 



Bm z , 



— m H m 
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The last term implements the constraint of charge conserva- 
tion with the Lagrange multiplier — qVo. From the Euler equa- 
tion dh/d(Sn) — we obtain 



on = J \ v ~ v o) + "ir~2 "TT m ■ ( 9 ) 
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The constraint Sn = leads to 



V — Vq = TO . 

2q 



(10) 



Here, Vo is a reference potential, which is irrelevant for the 
measured voltages. From dh/dm = together with Sn = 
we obtain m x — m y — and the standard result a'(T — 
T c )m z + f3m? z — B = 0. The elementary solution for both 
T < T, and T > T c is 
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where b = B/B s . For weakly inhomogeneous magnetic 
fields we obtain the linear response of the voltage to the mag- 
netic field gradient by expanding Eq. dl0> . 
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An analytical expression can be obtained from Eq. M \\ . Using 
dimensionless quantities, we find 
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with V s = 7/(roV/3 £ o e )- The sign is given by the sign of 
the carrier charge. Results for different magnetic fields are 
shown in Fig. ^ In me limit of vanishing magnetic field we 
obtain, in the ferromagnetic phase, dV/dB (V s /B s )~ l = ± 
(g/2) l/y/^-t. The linear-response coefficient thus diverges 
as T c is approached from below. 

We next estimate the order of magnitude of dV/dB. While 
the Landau theory does not apply down to T = 0, extrapolat- 
ing to T = results in a magnetization of the correct order of 
magnitude. For a rough estimate we thus write 
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dB 
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Here, rj w 5.4 x 10 25 Km 3 has been estimated 17 from Ref.H, 
for the magnetization we take the maximum value for 5 % Mn 
in (Ga,Mn)As, m (T = 0) w 5.13 x 10 4 A/m, and T c is set 
to 110 K. This gives a typical scale of dV/dB\B=o ~ 7.9 x 
10~ 4 V/T at low temperatures and the coefficient increases 
like (T c — T)~ 1/2 as the temperature is increased. This result 
suggests that the effect should be clearly measurable. 

In the paramagnetic phase, the magnetization to is, to 
leading order, linear in magnetic field. Thus the voltage, 
Eq. fllOi , is quadratic in B and there is no linear response 
at zero field. Right at the critical point the Landau-theory 
result for the magnetization is m z — (B//3) 1/3 , leading to 
V = (a , T 1 /2q)(B/P) 2 ^. 

Next, we briefly discuss the voltage due to a large change 
in magnetic field. It may be possible to measure this voltage 
by first grounding a metallic DMS sample in zero magnetic 
field, removing the ground connection, and then applying a 
strong, uniform field. Since higher magnetization attracts 
carriers, a negative (positive) potential difference between 
sample and ground is expected for positive (negative) car- 
rier charge. Equation dlOt gives for vanishing magnetic field 
— Vq = {a' if) I (2q) TOq, where too is the (uniform) magnetiza- 
tion in zero magnetic field, which is mo = \J ' —a'{T — T c ) / [3 
in the ferromagnetic phase and too = in the paramagnetic 
phase. Consequently, 
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V is the potential that would be necessary to maintain Sn = 
if the sample were still connected to a charge reservoir. If 
the sample is isolated, a voltage — V with the reverse sign is 
measured. In terms of dimensionless quantities we obtain 



V 
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where the upper (lower) sign applies to positive (negative) 
carrier charge. For small magnetic fields we get the same 
linear-response result as above. The full nonlinear response 
is obtained by inserting Eq. (1111 into Eq. d!5l >. The result- 
ing voltage is plotted as a function of magnetic field for var- 
ious temperatures in Fig. [2] In the ferromagnetic phase, the 
voltage crosses over from a linear B dependence to £? 2/3 at 
B/B s ~ < 3/2 , see Eq. (lilt . In the paramagnetic phase there 
is a crossover from B 2 to B 2 73 at the same scale. The inset in 
Fig-Elshows the dependence of \V\ on temperature at a finite 
magnetic field. The cusp at T = T c stems from the zero-field 
magnetization too in Eq. d!5l >. 

All results up to this point have been obtained from Landau 
mean-field theory. The question arises whether the results sur- 
vive if fluctuations are ta ken i nto account. After all, the aver- 
age square magnetization to 2 becomes an analytical function 
of temperature through T c . 

We restrict ourselves to temperatures T » T c and to small 
magnetic fields. Let us assume that charge fluctuations are 
fast compared to fluctuations in the magnitude of m. This 
seems reasonable since charge density fluctuations have a typ- 
ical timescale of the inverse plasma frequency. In this case the 




FIG. 2: (Color online) Magnitude of voltage —V between sample 
and ground as a function of magnetic field for various temperatures 
for a DMS sample in the uniform ferromagnetic or paramagnetic 
phases. Here, t = a'(T — T c )rQ/j is the dimensionless reduced 
temperature. If the voltage is scaled with 1/g, as it is plotted here, 
the curves do not depend on g. Inset: Voltage as a function of t for 
B/B s = 1. 



carrier concentration instantaneously follows the magnetiza- 
tion fluctuations. Equation dl 2i then becomes 



dV a'rj dm 2 
dB~^2q~ ~dB' 



(17) 



Equation Q shows that the averages of m z and to 2 are 
obtained from the exact free energy density / by toJ = 
—df/dB and to 2 = 2df/da. This implies the Maxwell 
relation 



dm 2 
~~dB 



da 



(18) 



For vanishing magnetic field and T < T c , one finds m7 oc 
(T c — T)— , where is the usual critical exponent (not to be 
confused with the coefficient (3 in the Landau functional) and 
T c is the true Curie temperature, taking fluctuations into ac- 
count. Since a is a linear function of T we obtain 



dV 
dB 



oc (T c - T)E 



0-i 



(19) 



B=0 



This quantity still diverges as T c is approached from below. 
Note that for the three-dimensional Heisenberg model, (3 ~ 
0.37. 21 Fluctuations thus make the divergence stronger. In the 
paramagnetic phase, toJ of course vanishes for B = and we 
still do not find a linear response. 



B. Stripe phase 

Below T* the Landau theory predicts an equilibrium 
state with periodic modulations of magnetization and carrier 
concentration. 17 We briefly review the pertinent properties of 



4 





1 1 






















-1 


1 


^1- 





mlm 




-1 1 

reduced temperature g 2 t 

FIG. 3: (Color online) Diagram showing the equilibrium phases as 
a function of reduced temperature g 2 t in a nonzero magnetic field 
B/B s = 1. For t > t* = a'(T* - T c )rf,/-y the equilibrium state 
is uniform with the magnetization given by the heavy solid line. For 
t < t* the magnetization shows periodic modulations spanning the 
cross-hatched magnetization interval. The lower limit of the mod- 
ulations is m = m s i n g. The heavy dashed line shows the magne- 
tization of the uniform solution for t < t* , which exists but has 
higher energy. The thin lines denote the uniform magnetization and 
upper limit of modulations for B — 0, for comparison. The left in- 
set shows the magnetization as a function of position for g 2 t = —2 
and B/B s — 1. The right inset shows the effective potential fto 
as a function of magnetization for the same parameter values. The 
double-headed arrow denotes the magnetization modulation. The de- 
tailed numerical values depend on the choice of g (g = 1 in this plot) 
and of B/B 3 , but the topology of the diagram does not. 



this state. The first two terms in the Landau functional H m , 
Eq. Q, can be interpreted as a potential for a fictitious particle 
with coordinates m. The uniform equilibrium states are deter- 
mined by the minima of this potential. In the stripe phase the 
coefficients are renormalized, 17 leading to the energy density 
h = (a/2) m 2 + 0/4) to 4 with 



/2 2 

a (1 — l c ) H „ r, 



2e 2 r 2 



/2 2 
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2e 2 r 2 ' 



(20) 
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where to 2 has to be calculated selfconsistently. It turns out 
that a modulation of to about the minimum of ho is stable 
and energetically favorable for m > m S i ng , where m S j ng = 
^e/a'rjy/eoe = m s /g. Modulations occur between two 
magnetizations with equal potential ho(m). For weak cou- 
pling, i.e., small rj or g, m s i ng is larger than the magnetization 
at the minimum of ho(m) and stable modulations are not pos- 
sible. For larger coupling, m S i ng is smaller than m at the min- 
imum, allowing modulations to occur. The minimum-energy 
solution is obtained for the maximum possible modulation am- 
plitude. In this case the lower turning point of the modulation 
approaches m s i ng »Ai 

In the present case we add the Zeeman term —Bm z to the 
energy density, which tilts the effective potential, as shown 
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FIG. 4: (Color online) Linear response dV/ dB for small magnetic 
field, as a function of reduced temperature. The jump in dV/dB at 
t — t* is discussed in the text. For t < t* , the result for the unstable 
uniform state would continue along the dashed line. The numerical 
calculation has been performed for the coupling g = 1. Inset: Phase 
diagram in the temperature-magnetic-field plane. The heavy solid 
line denotes the first-order transition of the ferromagnet, ending in a 
critical point at T c (i.e., t — 0). The thin solid lines emanating from 
t — t* (B = 0) denote the second-order transition between uniform 
and stripe phases. 



in the right inset in Fig. [3] Numerical integration shows that 
at any temperature the modulation with maximum amplitude 
still has the lowest energy. Figure [5] shows a typical phase 
diagram in the temperature-magnetization plane. A typical 
minimal-energy solution is shown in the left inset. 

The field-dependent transition temperature T*(B) to the 
stripe phase is determined by the condition m = m S j ng , where 
to is the uniform magnetization, Eq. dl 11 1. The solution is 



T* = T r 
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Using Eqs. @-(|6), this is equivalent to 

a'(T* - T c )r 2 1 
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(23) 



Thus the transition temperature depends linearly on B for ar- 
bitrary B. Of course, for large magnetic fields Landau the- 
ory becomes inapplicable. Note that while there is no sharp 
paramagnet-ferromagnet transition in nonzero field, the tran- 
sition to the stripe phase remains sharp. The reason is that 
translational symmetry is preserved by the applied field and 
thus can be spontaneously broken at this transition. The inset 
in Fig.f4]shows the resulting phase diagram. 
The analogue of Eq. dlOt in the stripe phase is 



V 
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If fluctuations are neglected, the average is a spatial average 
over the static modulation. The linear response dV/dB = 
(a'r)/2q) dm 2 /dB is calculated numerically and the result for 



B = is shown in Fig. |4] The linear-response coefficient 
shows a discontinuity at T* within Landau theory. This phe- 
nomenon is accompanied by a similar jump in the susceptibil- 
ity at T* (not shown). Physically, the system is more easily 
polarizable in the stripe phase, since it has an additional pa- 
rameter, the wavelength A, that can be adapted to the applied 
field. In agreement with this interpretation, the inset in Fig. [4] 
shows that the stripe phase is stabilized by a magnetic field. 



III. THERMOPOWER 

Since the magnetization changes with temperature and is 
coupled to the carrier concentration, there should also be a 
change in carrier concentration or potential with temperature. 
In the case of an electrically isolated sample one expects a 
nonzero thermopower 



Q 



dV 
df 



(25) 
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Consider the following setup: One end of a DMS sample is 
kept at a temperature T\ , the other at temperature T2 > T\ . 
A discussion of the standard origin of the thermopower can 
be found in textbooks. 22 The result is that the Q has the same 
sign as the carrier charge q. For this reason the thermopower 
is often meassured to obtain the sign of this charge in semi- 
conductors, including DMS. '^12 

Here, we predict an additional contribution in DMS: In the 
same setup, with T\ < T2 < T c , the magnetization m is larger 
in the cool region. Carriers are attracted to high-magnetization 
regions, leading to a current which results in the accumulation 
of positive (negative) charge in the cool region for positively 
(negatively) charged carriers. This charge generates an elec- 
tric field that prevents further charge flow. The electric poten- 
tial <p is higher (lower) in the cool region for positive (nega- 
tive) carrier charge q. With the explicit minus sign in Eq. (125 \ 
this leads to a thermopower of the same sign as the normal 
contribution. Proceeding as in Sec.||T]we obtain 



ct'rj 



V = 4>(T 2 ) - <j»{Ti) = m 2 (T 2 )-m 2 (Ti 
Zq 
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We restrict ourselves to vanishing magnetic field. In the uni- 
form phase we can write, neglecting fluctuations, to 2 = m 2 . 
In Landau theory, Eq. (I26> gives 



Q = Q V 



/2 

a r/ 
2qJ' 



(27) 



Thus the thermopower is independent of temperature for 
T* < T < T c . On the other hand, for T > T c there is no 
magnetization and this contribution to the thermopower van- 
ishes. We thus find a jump at T c . 

In the stripe phase, for T < T*, we find numerically that 
the thermopower increases in magnitude with temperature. At 
T* it shows a downward jump, since the average square mag- 
netization to 2 changes slope as a function of temperature at 
the transition. Thermopower and m 2 are shown in Fig.|5] 
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FIG. 5: (Color online) Left axis: Thermopower Q in units of the con- 
stant thermopower Q U ni in the uniform phase for vanishing magnetic 
field, as a function of reduced temperature. Also shown is the aver- 
age square magnetization m' 2 (right axis). The numerical calculation 
for t < t* has been performed for g = 1. 



To estimate the order of magnitude of the thermopower, 
we again replace a'//3 by m 2 (T = Q)/T c , leading to Q = 
—(3r]m A (T = Q)/2qT 2 . The only additional parameter needed 
beyond those given in Sec.|II|is (3. By comparing the mean- 
field result for the gain in energy density due to magnetic or- 
dering of the Heisenberg model on a simple cubic lattice to 
the corresponding gain for the Landau theory, we obtain 



/3 



1 



6k B T c 



(,?/xb) 4 S 3 (S- 



2.1 x 1(T 12 Jm/A 4 , (28) 



where g w 2 is the g-factor for the impurity spins with S = 
5/2 and riMn is the concentration of impurities, again assumed 
to be 5 % of cation sites. This gives Q w 2.0 x 1CT 3 V/K. 
In nonmagnetic semiconductors showing thermally activated 
conduction, the magnitude of the thermopower is of the order 
of 1(T 5 to 1CT 4 V/K. In metals it is a few times 1CT 6 V/K. 
Thus the thermopower generated by the carrier-magnetization 
coupling is expected to be larger than the normal contribution. 

The thermopower has been measured for various DMSp&A^ 
but usually only to infer the sign, which is not changed by 
the physics discussed here. A study of the thermopower in 
(Ga,Mn)As is under way. 20 It would be worthwhile to look 
for an anomalously large thermopower in DMS. However, the 
main signature of the magnetization-induced effect would be 
a downward jump at T c . 

Finally, we briefly turn to the qualitative effect of fluctua- 
tions: The average square magnetization becomes analytic in 
temperature through T c so that we expect the jump in Q to be 
replaced by a continuous crossover. Its width should be given 
by the Ginzburg criterion for the fluctuation-dominated tem- 
perature interval. In real DMS, disorder plays an important 
role|S which may lead to additional broadening. However, for 
high-quality metallic samples of (Ga,Mn)As, the magnetiza- 
tion curves close to T c show a sharp decrease^i&SiiS similar 
to the Heisenberg model on a regular lattice, suggesting that 
disorder is not dominant. 
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IV. SUMMARY AND CONCLUSIONS 

The strong coupling between carrier concentration and 
magnetization in ferromagnetic DMS has been found to lead 
to a dependence of the electric potential on the magnetization. 
If the magnetization is controlled by an applied magnetic field, 
a magnetic-field-induced voltage is expected. Two possible 
setups are discussed in this paper, which allow to study the 
linear response for weakly nonuniform magnetic fields and 
the nonlinear voltage induced by a strong field. The linear- 
response coefficient in the limit of small field is estimated to 
be of the order of 10~ 3 V/T for (Ga,Mn)As at low tempera- 
tures. It diverges as T c is approached from below. This singu- 
larity persists if fluctuations are taken into account. 

On the other hand, if the magnetization is changed by vary- 
ing the temperature, the same physics leads to a variation 
of electric potential with temperature, i.e., a thermopower. 
Within Landau theory, the thermopower is temperature- 
independent in the uniform ferromagnetic phase and shows 
a discontinuity at the Curie temperature T c , which is smeared 
out by fluctuations. For (Ga,Mn)As the estimate for this con- 



tribution to the termopower is of the order of 10 _3 V/K, 
which is what one would have guessed from the magnitude 
of the magnetic-field-induced voltage. All these results are 
obtained for a uniform equilibrium magnetization. 

A nonuniform, stripe-like equilibrium magnetization is 
possible in DMS at lower temperatures, T < J 1 *. 17 The phase 
diagram of uniform and stripe phases is obtained within Lan- 
dau theory. If the stripe phase were realized, the magnetic- 
field-induced voltage would be qualitatively similar to the uni- 
form case, but the linear response dV/ dB would show a dis- 
continuity at T*. The effect is larger in the stripe phase be- 
cause the magnetization can adapt more easily to the applied 
magnetic field by varying the wavelength of the stripe pat- 
tern. The thermopower in the stripe phase is temperature- 
dependent, unlike in the uniform phase, and shows a down- 
ward jump at T*. 
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